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CONFORMAL MODULUS ON DOMAINS WITH STRONG 
SINGULARITIES AND GUSPS 

KARRI HAKULA*, ANTTI RASILA^, AND MATTI VUORINEN* 

Abstract. We study the problem of computing the conformal modulus of rings and quadrilater¬ 
als with strong singularities and cusps on their boundary. We reduce this problem to the numerical 
solution of the associated Dirichlet and Dirichlet-Neumann type boundary values problems for the 
Laplace equation. Several experimental results, with error estimates, are reported. In particular, 
we consider domains with dendrite like boundaries, in such cases where an analytic formula for the 
conformal modulus can be derived. Our numerical method makes use of an hp-FEM algorithm, 
written for this very complicated geometry with strong singularities. 
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1. Introduction. The conformal modulus is an important tool in geometric 
function theory [1], and it is closely related to certain physical quantities which also 
occur in engineering applications. For example, the conformal modulus plays an 
important role in determining resistance values of integrated circuit networks (see 
e.g. [MIIID]). We consider both simply and doubly-connected bounded domains. By 
definition such a domain can be mapped conformally either onto a rectangle or onto 
an annulus, respectively. For the numerical study of these two cases we dehne the 
modulus h as follows. In the simply connected case, we hx four points on the boundary 
of the domain, call a domain with these hxed boundary points a quadrilateral, and 
require that these four points are mapped onto vertices (0,0), (1,0), (l,h), (0,h) of 
the rectangle. In the doubly connected case we require that the annulus is {{x,y) : 
exp(—h) < + y'^ < 1} . Doubly connected domains are also called ring domains or 

simply rings. Surveys of the state of the art methodologies in the field are presented 
in the recent books by N. Papamichael and N. Stylianopoulos [36] and by T. Driscoll 
and L.N. Trefethen [ISj. Various applications are described in nanHiiiniiiT]. In the 
past few years quadrilaterals and ring domains of increasing complexity have been 
studied by several authors 

We consider the problem of numerically determining the conformal modulus on 
certain ring domains with elaborate boundary. Due to the structure of the bound¬ 
ary, the problem is numerically challenging. On the other hand, the ring domain is 
characterized by a triplet of parameters (r, to,p), its construction is recursive, and 
yet its conformal modulus can be explicitly given. Varying the parameter values or 
the recursion level of the construction one can increase the computational challenge 
and therefore this family of domains forms a good set of test problems. In particular, 
error estimates can be given. For a figure of a domain in this family see Figure |3.1| 
The boundaries of these ring domains are point sets of dendrite type (i.e. continua 
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without loops) [48]. We apply here the /ip-FEM method developed in [22] for the 
computation of the moduli of these ring domains and report the accuracy of our al¬ 
gorithm. Furthermore, we use the algorithm of |21] to numerically approximate the 
canonical conformal mapping of the above class of domains (see Figure 3.2 (c)). The 
conjugate problem of the original ring problem solved in the approximation of the 
conformal mapping can be interpreted as a simplified crack problem. 


Due to the pioneering work of I. Babuska and his coauthors [Sll6l[7] the opti¬ 
mal convergence rate of the /ip-FFM method is well studied and experimentally also 
demonstrated in some fundamental basic situations. Our main results are the compu¬ 
tation the moduli of ring domains, given in the form of various error estimates, where 
we compare three error norms in the case of rings with elaborate boundary: (1) Fxact 
error (2) Auxiliary space estimate, based on /ip-space theory (3) So called reciprocal 
error estimate, introduced in |22j . We also attain the nearly optimal convergence in 
accordance with the theory of [Slll[7]. All these three error estimates behave in the 
same way. Because this is the case, there is some justification to use the estimates 
(2) and (3) also in the common case of applications when the exact value of the mod¬ 
ulus is not know and hence estimate (1) is not available. We apply our methods to 
compute the modulus of a quadrilateral considered by Bergweiler and Eremenko [10] . 
The boundary of this domain has cusp-like singularities. 


1.1. Project background and history. This paper is a culmination of a ten- 
year research project, arising from questions related to the work of Betsakos, Samuels- 
son and Vuorinen m- The original goal of the project was to develop accurate nu¬ 
merical tools suitable for studying effects of geometric transformations in function 
theory (see e.g. [El [24]). Experimental work towards this goal was carried out by 
Rasila and Vuorinen in the two small papers [37l EH] , and further work along the same 
lines was envisioned. However, it was quickly discovered that the AFEM package of 
Samuelsson used in the above papers is not optimal for studying very complex ge¬ 
ometries arising from certain theoretical considerations, as the number of elements 
used in such computations tends to become prohibitively large. This problem led us 
into the higher order /ip-FEM algorithm implemented by Hakula. Efficiency of this 
method for numerical computation of conformal modulus had been established in the 
papers [HH] and [53], the latter of which deals with unbounded domains. Recently, an 
implementation of this algorithms for the purpose of numerical conformal mapping 
was presented in [5T] . 

2. Preliminaries. In this section central concepts to our discussion are intro¬ 
duced. The quantities of interest from function theory are related to numerical meth¬ 
ods, and the error estimators arising from the basic principles are defined. 

2.1. Conformal Modulus. A simply-connected domain D in the complex plane 
C whose boundary is homeomorphic to the unit circle, is called a Jordan domain. 
A Jordan domain D, together with four distinct points zi, Z 2 , Z 3 , in dD, which 
occur in this order when traversing the boundary in the positive direction, is called a 
quadrilateral and denoted by {Di; zi, Z 2 , Z 3 , 04 ). If /: D ^ fD is a conformal mapping 
onto a Jordan domain fD, then / has a homeomorphic extension to the closure D 
(also denoted by /). We say that the conformal modulus of {D; zi, Z 2 , Z 3 , Z 4 ) is equal 
to h > 0, if there exists a conformal mapping f ot D onto the rectangle [0,1] x [0, h], 
with f{zi) = l + ih, /(Z 2 ) = ih, /(zs) = 0 and /(Z 4 ) = 1. 

It follows immediately from the definition that the conformal modulus is invariant 
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under conformal mappings, i.e., 

zi,Z 2 ,Z 3 , Z 4 ) = M(/D; f(zi), /(Z 2 ), /(zs), f(z 4 )), 


for any conformal mapping /: D f{D) such that D and f{D) are Jordan domains. 

For a curve family F in the plane, we use the notation M(F) for its modulus 
[30] , For instance, if F is the family of all curves joining the opposite 6-sides within 
the rectangle [0, a] x [0,6], a, 6 > 0, then M(r) = 6/a. If we consider the rectangle 
as a quadrilateral Q with distinguished points a + ib,ib,0,a we also have M{Q-,a + 
ib, ib, 0,a) = b/a , see [T1I3D]. Given three sets D, E, F we use the notation A(E, F; D) 
for the family of all curves joining E with F'm D . 

2.2. Modulus of a quadrilateral and Dirichlet integrals. One can ex¬ 
press the modulus of a quadrilateral {D; zi, Z2, z^^ Z4) in terms of the solution of the 
Dirichlet-Neumann problem as follows. Let 7 j, j = 1, 2, 3,4 be the arcs of dD between 
(z 4 , zi) , {zi, Z2) , {z2, Z3) , ( 2 : 3 , Z4), respectively. If u is the (unique) harmonic solution 
of the Dirichlet-Neumann problem with boundary values of u equal to 0 on 72 , equal 
to 1 on 74 and with dujdn = 0 on 71 U 73 , then by [TJ p. 65/Thm 4.5]: 


M{D;zi,Z 2 ,Z 3 ,Z 4 ) = // \\7u\'^dxdy 


( 2 . 1 ) 


D 


The function u satisfying the above boundary conditions is called the potential func¬ 
tion of the quadrilateral (D; zi, 22 , 23 , 24 ). 

2.3. Modulus of a ring domain and Dirichlet integrals. Let E and F be 

two disjoint compact sets in the extended complex plane Coo. Then one of the sets 
E, F is bounded and without loss of generality we may assume that it is i?. If both 
E and F are connected and the set R = Coo \{EU F) is connected, then R is called 
a ring domain. In this case i? is a doubly connected plane domain. The capacity of 
R is defined by 



where the infimum is taken over all nonnegative, piecewise differentiable functions 
u with compact support in i? U FI such that u = 1 on E. It is well-known that 
there exists a unique harmonic function on R with boundary values 1 on E and 0 
on F. This function is called the potential function of the ring domain i?, and it 
minimizes the above integral. In other words, the minimizer may be found by solving 
the Dirichlet problem for the Laplace equation in R with boundary values 1 on the 
bounded boundary component E and 0 on the other boundary component F. A ring 
domain R can be mapped conformally onto the annulus {2 : < jzj < 1}, where 

M — M(i?) is the conformal modulus of the ring domain R. The modulus and capacity 
of a ring domain are connected by the simple identity M(i?) = 27r/capi?. For more 
information on the modulus of a ring domain and its applications in complex analysis 
the reader is referred to [nnsiiMii^- 

2 . 4 . Hyperbolic Metrics. The hyperbolic geometry in the unit disk is a pow¬ 
erful tool of classical complex analysis. We shall now briefly review some of the main 
features of this geometry, necessary for what follows. First of all, the hyperbolic 
distance between a:, y G D is given by 
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In addition to the unit disk D, one usually also studies the upper half plane IHI as a 
model of the hyperbolic geometry. For a;, y € El we have {x = {xi,X 2 )) 

pmix^y) = arcosh 1 + f-- 

V 2x2y2 

If there is no danger of confusion, we denote both pa(z,w) and poiz,w) simply by 
p{z,w). We assume that the reader is familiar with some basic facts about these ge¬ 
ometries: geodesics, hyperbolic length minimizing curves, are circular arcs orthogonal 
to the boundary in each case. 

Let zi, Z 2 , Z 3 , Z 4 be distinct points in C. We define the absolute (cross) ratio by 


\ZI,Z2,Z3,Z4\ 


Izi - 2 : 3 ! IZ 2 - Z4I 
l-Zl - Z2l |Z3 - Z4l ' 


( 2 . 2 ) 


This definition can be extended for zi, Z 2 , Z 3 , Z 4 € Coo by taking the limit. An im¬ 
portant property of Mobius transformations is that they preserve the absolute ratios, 

i.e. 


If(zi),f(z2),f(z3),f(z4)l = |zi,Z2,Z3,2^4|, 

if /: Coo —>■ Coo is a Mobius transformation. In fact, a mapping /: Coo —>■ Coo is a 
Mobius transformation if and only if / it is sense-preserving and preserves all absolute 
ratios. 

Both for (D, pju) and (H, pjj) one can define the hyperbolic distance in terms of the 
absolute ratio. Since the absolute ratio is invariant under Mobius transformations, 
the hyperbolic metric also remains invariant under these transformations. In partic¬ 
ular, any Mobius transformation of D onto El preserves the hyperbolic distances. A 
standard reference on hyperbolic metrics is [2]. 

2.5. hp-FEM. In this work the natural quantity of interest is always related 
to the Dirichlet energy. Of course, the finite element method (FEM) is an energy 
minimizing method and therefore an obvious choice. The continuous Galerkin hp- 
FEM algorithm used throughout this paper is based on our earlier work [52]. Brief 
outline of the relevant features used in numerical examples below is: Babuska-Szabo - 
type p-elements, curved elements with blending-function mapping for exact geometry, 
rule-based meshing for geometrically graded meshes, and in the case of isotropic p 
distribution, hierarchical solution for all p. The main new feature considered here is 
the introduction of auxiliary subspace techniques for error estimation. 

For the types of problems considered here, theoretically optimal conforming hp- 
adaptivity is hard. The main difficulty lies in mesh adaptation since the desired 
geometric or exponential grading is not supported by standard data structures such 
as Delaunay triangulations. Thus, the approach advocated here is a hybrid one, where 
the problem is first solved using an a priori /ip-algorithm after which the quality of 
the solution is estimated using error estimators specific both for the problem and 
the method, provided the latter are available. For instance, the exact solution or for 
problems concerning the conformal modulus the so-called reciprocal error estimator. 
The a priori algorithm is modified if the error indicators suggest modifications. If this 
occurs, the solution process is started anew. 

In the numerical examples below the computed results are measured with both 
kinds of error estimators giving us high confidence in the validity of the results and 
the chosen methodology. 
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2.5.1. Auxiliary Subspace Techniques. Consider the abstract problem set¬ 
ting with Vas the standard piecewise polynomial hnite element space on some dis¬ 
cretization T of the computational domain D. Assuming that the exact solution 
u G Hq{D) has finite energy, we arrive at the approximation problem: Find u G V 
such that 


a{u,v) = l{v) (= a{u,v)), G V, (2-3) 

where «(■,•) and /(•), are the bilinear form and the load potential, respectively. Ad¬ 
ditional degrees of freedom can be introduced by enriching the space V. This is 
accomplished via introduction of an auxiliary subspace or “error space” W C Hq (D) 
such that V nW = {0}. We can then define the error problem: Find e G V such that 

a(e, v) = l{v) — a{u, v){= a(u — u, v)), Vw G W. (2.4) 


In 2D the space W, that is, the additional unknowns, can be associated with element 
edges and interiors. Thus, for hp-methods this kind of error estimation is natural. 
For more details on optimal selection of auxiliary spaces, see [20] . 


The solution e of (2.4) is called the error function. It has many useful properties 


for both theoretical and practical considerations. In particular, the error function can 
be numerically evaluated and analysed for any finite element solution. This property 
will be used in the following. By construction, the error function is identically zero 
at the mesh points. In Figure [373] one instance of a contour plot of the error function 
(with a detail) is shown. This gives an excellent way to get a qualitative view of the 
solution which can be used to refine the discretization in the hp-sense. 

Let us denote the error indicator by a pair (e, &), where e and b refer to added 
polynomial degrees on edges and element interiors, respectively. It is important to 
notice that the estimator requires a solution of a linear system. Assuming that the 
enrichment is fixed over the set of p problems, it is clear that the error indicator 
is expensive for small values of p but becomes asymptotically less expensive as the 
value of p increases. Following the recommendation of m, our choice in the sequel 
is (e, 6) = (1,2) unless specified otherwise. 

Remark 2.1. In the case of (Q,b)-type or pure bubble indicators, the system 
is not connected and the elemental error indicators can be computed independently, 
and thus in parallel. Therefore in practical cases one is always interested in relative 
performance of {0,b)-type indicators. 

2.6. Reciprocal Identity and Error Estimation. Let Q be a quadrilateral 
defined by points zi,Z 2 ,Z 3 ,Z 4 and boundary curves as in Section [2.1| above. The 
following reciprocal identity holds: 


As 


M{Q-,zi,Z2,Z3,Z4)M{Q-,Z2,Z3,Zi,zi) = 1 . 

we shall use the test functional 

Zi, Z2, Zs„ Z 4 )U{Q-, Z2, Z 3 ., Zi, Zi) - l| 


(2.5) 


( 2 . 6 ) 


which by (2.5) vanishes identically, as an error estimate. 

As noted above, the error function e can be analysed in the sense of FEM- 
solutions. Our goal is to relate the error function given by auxiliary space techniques 
and the reciprocal identity arising naturally from the geometry of the problem. Let 
us first define the energy of the error function e as 


f(e) = JJ \Ve\^dxdy. 


(2.7) 
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Using (2.71 the reciprocal error estimation and the error function e introduced 
above can be connected as follows: Let Oi and 02 be the moduli of the original and 
conjugate problems, ei, £ 2 , and ei = f (ei), £2 = ^(£ 2 ), the errors and their energies, 
respectively. Taking e = max{|ei|, |e 2 |} we get via direct computation: 


|1 — (oi + ei)(a 2 + £ 2)1 ^ |aiC 2 + 0 , 2^1 + €i£ 2 | ^ 2 £max{ai, l/oi} + O(i^). (2.8) 


Neglecting the higher order term one can solve for i and compare this with the esti¬ 
mates given by the individual error functions. 

3. The Dendrite. A compact connected set in the plane is called dendrite-like if 
it contains no loops. We introduce a new parametrized family of ring domains whose 
boundaries have dendrite-like boundary structure and whose modulus is explicitly 
known in terms of parameters. In numerical conformal mapping one usually considers 
domains whose boundaries consist of finitely many piecewise smooth curves. Very 
recently in |39j these authors considered conformal mapping onto domains whose 
boundaries have “infinitely many sides”, i.e., are obtained as result of a recursive 
construction. One of the examples considered in |39| was the domain whose boundary 
was the von Koch snowflake curve. 

In this section we will give a construction of a ring domain whose complementary 
components are C \ D and a compact connected subset C (r, p, m) of the unit disk D 
depending on two positive integer parameters p, m and a real number r > 0. The set 
C{r,p,m) consists of finitely many pieces, each of which is a smooth curve, and the 
set is acyclic, i.e., does not contain any loops. The number of pieces is controlled by 
the integers (p, m) and can be arbitrarily large when p and m increase. 

3.1. Theory. Recall that the Grotzsch ring Rcir) = D \ [0,r], r G (0,1) has the 
capacity cap(i?G(?')) = 27r/p(r), were p(r) is the Grotzsch modulus function (cf. [21 
Ghapter 5]): 




dx 


-^(1 — a;2)(l — r^a;^) 


with usual notation r' = y/l — r^. Let r G (0,1), and let Dr =Bi \ ([—r, r]U [—ir, ir]). 
The conformal mapping f(z) = ^—z maps the Grotzsch ring Rcif) (excluding the 
positive) real axis onto the sector {z : \ argzj < 7r/2}. Let Ur be the potential function 
associated with Rq^x). Then, by symmetry, it follows that the potential function 
Ur o f can be extended to the domain Dr by Schwarz symmetries so that it solves 
the Dirichlet problem associated with the conformal capacity of Dr- It follows that 
cap{Dr) = 8Tr/p{r‘^). Obviously, the similar construction is possible for any integer 
TO > 3. 

One may continue the process to obtain further generalizations. Start with a 
generalized Grotzsch ring with to > 3 branches. Choose one of the vertices of the 
interior component. Map this point to the origin by a Mobius automorphism of the 
unit disk. Make a branch of degree p > 2 to the origin by using the mapping z 1 —^ z^^^, 
and extend the potential function to the whole disk by using Schwarz symmetries. The 
resulting ring has capacity 27rTOp/p(r'"). An example of the construction is given in 
Figure [XT] 

Again, it is possible to further iterate the above construction to obtain ring do¬ 
mains with arbitrarily complex dendrite-like boundaries. Let to > 3, M > 1, and 
let be integers such that pj > 2 for all j = 1,2,..., M. For each j = 1,2,..., M, 
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choose one of the vertices Zj of the interior component. Let Wj be the point on the 
line tzj, t > 0 so that \wj\ = 1. We may assume that Zj < 0, Wj = —1 and the line 
segment [— 1 , 2 ^] does not intersect with the interior component except at the point 
Zj. Map the point Zj to the origin by Mobius automorphism gj of the unit disk so 
that gj{—\) = —1. Now map the domain D \ [—1, 0] onto the symmetric disk sector 
by the mapping hj{z) = , and extend the potential function to the whole disk 

by using Schwarz symmetries. By repeating this construction for all j = 1,2,..., M, 
we obtain a ring domain with the conformal capacity 


27rmpip2 ■■■PM 


3.2. Numerical Experiments. We consider two cases described in Table 3.1 
and Figures and Using the /ip-refinement strategy at the tips of the dendrite 
and the inner angles (120°) we obtain exponential convergence in the reciprocal error 


(Figure 3.4). In Figure 3.3 we also show the error function of type (1,2) over the 


whole domain as well as a detail which clearly shows the non-locality of the error 
function. As expected the errors are concentrated at the singularities and in the 
elements connecting the singularities to the boundary. One should bear in mind, 
however, that the reciprocal errors are small already at p = 10 used in the figures. 

The effect of error balancing is 


The error estimates are shown in Figure 3.5 


evident. For the larger capacity, the reciprocal error conincides with the true error, but 
overestimates the smaller one. However, in both cases the rates are correct, only the 
constant is overly pessimistic. The auxiliary space error estimate underestimates the 
error slightly, again with the correct rate. This is exactly what we would expect, since 
increasing the polynomial orders in the scheme should increase the error estimate, if 
the underlying solution has converged to the correct solution. 


Table 3.1: Tests on dendrite-problems. The errors are given as Klog^o |error|]|. 


Case 

Parameters 

Method 

Errors 

Sizes 

M{Qs) 

1 

r = 1/20, m = A, p = 3 

hp, p = 16 

9 (9) 

159865 (160161) 

5.63968609980242 

2 

r = 1/20, m = 5, p = 7 

hp, p = 12 

9 (9) 

199921 (200809) 

13.437951766839522 


4. Domains with Cusps. 

4.1. Hyperbolic Quadrilateral. Let Qs be the quadrilateral whose sides are 
circular arcs perpendicular to the unit circle with vertices e*®, and 

e“®®. We call quadrilaterals of this type hypebolic quadrilaterals as their sides are 
geodesics in the hyperbolic geometry of the unit disk. We approximate values of the 
modulus of Qs- 

Next we determine a lower bound for the modulus of a hyperbolic quadrilateral. 
Let 0<a</3<7<27r. The four points 1, e®“, e®^, e®^ determine a hyperbolic 
quadrilateral, whose vertices these points are and whose sides are orthogonal arcs 
terminating at these points [301 HI]. We consider the problem of finding the modulus 
(or a lower bound for it) of the family F of curves within the quadrilateral joining the 
opposite orthogonal arcs (e®“, e®^) and (e®^, 1) within the quadrilateral [30] . It is easy 
to see that we can find a Mobius transformation /i of D onto El such that h{l) = 1, 
/i(e®“) = t, ft,(e®^) = —t, /i(e®'’') = —I for some t > 1. The number t can be found 
by setting the absolute ratios |1, e®“, e®^, e®^| and |l,f, —t, I| equal, and solving the 
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(a) Generalized Grotzsch ring Ra{r,m): (b) Map the chosen point to the origin by a 
r = 1/4, m = 6. Mobius automorphism of the unit disk. 



(c) Make a branch of degree p = 7. (d) Extend the potential function to the 

whole disk. 


Fig. 3.1: Dendrite Construction: r = 1/4, m = 6, p = 7. 


resulting quadratic equation for t because Mobius transformations preserve absolute 
ratios. The image quadrilateral has four semicircles as its sides, the diameters of these 
are [—1,1], [1, t], [—t, t], [—t, —1] and the family h{r) has a subfamily A consisting of 
radial segments 

(j) e {e,Tr - e), sin0=|-^. 

Obviously, for 0 = 0 we obtain an upper bound. Therefore 

'jr 'TT _ Of) 

- -> M(/t(r)) > M(A) = --. 

logt - V V G - V 1 
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(c) Conformal mapping. 


Fig. 3.2: Dendrite 1: r = 1/20, m = 4, p = 3. 



Fig. 3.3: Dendrite 2: r = 1/20, m = 5, p = 7. 



Fig. 3.4: Dendrite 1: Reciprocal error; log-plot: Error vs p. 
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(a) Estimated error. (b) Estimated error (Conjugate). 

Fig. 3.5: Dendrite 1: Estimated errors; log-plot: Error vs p; Solid line = Reciprocal 
estimate, Dashed line = Auxiliary space estimate, Dotted line = Exact error. 



(a) Case 1: s = 7r/4. (b) Case 2: s = 3n/8. 

Fig. 4.1: Hyperbolic rectangles: Refined meshes. 


4.1.1. Numerical Experiments. Similarly as before, the examples of this sec¬ 
tion are outlined in Table [4H] and Figures |4T] |4. 2 1 The meshes are refined in exactly 
same fashion so that any differences in convergence stem only from the difference 
in the geometric scaling. As shown in Figure |4.3| the convergence in the reciprocal 
error is exponential, but with a better rate for the symmetric case. Moreover, for the 
symmetric domain both error estimates coincide. 

4.2. Other Quadrilaterals with Cusps. Next we consider quadrilaterals Q = 
Q{D; zi, Z 27 -^ 3 ; Z 4 ) where neither of the components of Cao\dD is bounded. The modu¬ 
lus of the following quadrilateral has been obtained by W. Bergweiler and A. Eremenko 
m, who studied this question in connection to an extremal problem of geometric 
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jugate). 

Fig. 4.2: Hyperbolic rectangles: Potential functions. 


10-2 



10-®-V- 

10 -®-- 

10 -^-X- 

10 -s-'X- 

10 -®-^ V- 

10-10 -- 

10-11 -X- 

1 0-12 I------- 

123456789 

(a) Case 1. 



(b) Case 2. 


Fig. 4.3: Hyperbolic rectangles: Reciprocal errors; log-plot: Error vs p. 


function theory introduced by A.A. Goldberg in 1973. 

4.2.1. Example I. Consider the strip the closed unit disk is removed: 

Hi = {z : —3 < Re 0 < 1} \ D. 

Let the four vertices Zj, j = 1, 2, 3,4 on the boundary of D be 1, oo, oo, 1, in counter¬ 
clockwise order. Then, all the angles at vertices are equal to 0. 

First we map the domain in question to a bounded domain so that the line 
{z : Rez = 1} maps to unit circle, and the real axis remains fixed. After the Mobius 
transformation we may assume that we are computing in the unit disk D. For con¬ 
venience, consider disks D(1 — t,t) and D(—1 -|- s,s) internally tangent to the unit 
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(b) Case 2. (c) Case 2 (Conjugate). 


Fig. 4.4: Hyperbolic rectangles: Estimated errors; log-plot: Error vs p; Solid line = 
Reciprocal estimate, Dashed line = Auxiliary space estimate. 


Table 4.1: Tests on hyperbolic quadrilaterals. The errors are given as | [logj^Q |error|] |. 


Case 

Method 

Errors 

Sizes 

M{Qs) 

1 

hp, p = 12 

12 

10225 

1 

2 

hp, p = IQ 

11 

17985 

3.037469188986459 


circle at the points 1 and 1, respectively, with s,t G (0,1/3). The corner points of the 
quadrilateral are 1, —1, —1,1, with a zero angle at each of the corners. We denote the 
respective radii of the disks by s and t. 

We have computed numerically the modulus of the family of curves joing the two 
disks within the domain £>2 = D \ (D(l — t, t) UD(—1 -|-s, s)). It is the reciprocal of the 
modulus of the family of curves joing the upper semicircle with the lower semicircle 
within the same domain. The results are summarized in Table 

An estimate for the case s = t = \/2 — 1 « 0.41421 is obtained by Bergweiler 
and Eremenko, with numerical values that agree with our results up to 6 significant 
digits, in m- The conformal modulus in this case is approximately 2.78234 (AFEM), 
or 2.7823418086 (hp-FEM, with error number = 10, p = 21). Our results agree with 
the result of m- 

4.2.2. Example II. Consider the domain D (a hexagon) in the upper half-plane 
obtained from the half-strip 

{z = a:-|-iy:0<x<l, 0< y}, 

by removing two half-disks 

Cl = 0(7/24,1/24), C 2 = 0(5/12,1/12), 

where 0(z, r) denotes the disk centered at z G C with radius r > 0. Note that CiOK = 
[1/4,1/3] and (72 O M = [1/3,1/2]. We compute the moduli of two quadrilaterals: 

Qi = {D- 00 ,0,1/2,1), Q 2 = {D- 0,1/4,1/2,1). 

Again, we hrst use the Mobius transformation 

2z - 1 


Z 


2z-b 1 
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Fig. 4.5: Quadrilaterals with Cusps: p-type Meshes. 


to map the domain D in question to a bounded domain. Then the boundary points 
of Qi map onto the points 1,—1,0,1/3, respectively. For <52, the boundary points 
are mapped onto the points —1,—1/3, 0,1/3. Quadrilaterals Qi and Q 2 , after the 
Mobius transformation, and the corresponding potential functions are illustrated in 
Figure [T6| 


4.2.3. Numerical Experiments. The numerical experiments differ from the 
previous cases since only the p-version is used. In other words, the meshes of Figure [475| 
are used as is, without any /i-refinement. As is evident in the convergence and error 
estimation graphs of Figures 4.7||4.9 n the cases where the local angles close to 7r/2, 
exponential convergence is achieved, but in the general case, when small geometric 
features are present, the convergence rates stall to algebraic and not exponential. 


Table 4.2: Tests on quadrilaterals with cusps. 


Case 

Method 

Errors 

Sizes 

M{Qs) 

1 

p, p = 16 

9 

1089 

2.7823418091539533 

2 

p, p = 16 

9 

1633 

1.8247899464782131 

3 

p, p = 16 

7 

2945 

1.7864319361374579 

4 

p, p = 16 

8 

2945 

0.8852475766134157 
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(a) Case 1. 



(b) Case 2. 




Fig. 4.6: Quadrilaterals with Cusps: Potentials. 


5. Conclusions. We have introduced a new class of ring domains, character¬ 
ized by three parameters, and given a formula for its modulus. By modifying the 
parameters, we obtain domains of increasing computational challenge. For some spe¬ 
cific sets of parameters we compute numerically the modulus and compare the true 
error and two error estimates and show that the these two error estimates behave in 
the same way as the true error. The computation is based on the /ip-FEM method 
and we show that nearly optimal convergence is obtained, when compared to the 
theory of I. Babushka and his coauthors. This class of domains could be used for 
benchmarking the numerical performance of FEM-software because of the scalability 
of computational challenge and the exactly known solution. 
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Fig. 4.7: Quadrilaterals with Cusps: Case 1: Reciprocal error, log-plot: Error vs p; 
Estimated error, log-plot: Error vs p; Solid line = Reciprocal estimate, Dashed line 
= Auxiliary space estimate. 
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